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Abstract 

We elucidate the role of time scales that determine heavy quark dynamics in the QGP. Quark- 
antiquark potentials extracted from the lattice are used to demonstrate the existence of open 
heavy flavor bound-state solutions in the vicinity of the critical temperature, and their light cone 
wavefunctions are obtained. We use these wavefunctions to calculate the in-medium modification 
of the heavy quark distribution and decay probabilities. For the case of high pj D or B mesons 
traversing the QGP, we combine the new meson formation and dissociation mechanism with the 
traditional parton-level charm and beauty quark quenching to obtain predictions for the heavy 
meson and non-photonic electron suppression in Cu+Cu and Pb+Pb collisions at RHIC and the 
LHC, respectively. 

1. Inadequacy of partonic level energy loss for heavy flavor suppression 

It is well known that partonic level energy loss of c and b quarks in the QGP is insufficient to 
explain the large suppression (Raa) of non-photonic e + + e~ measured at RHIC fl]]. To correctly 
assess the discrepancy between the expected heavy quark quenching and the measured Raa EL 
improved treatment of cold nuclear matter effects |13|,|4[ is necessary. It is, therefore, important to 
establish the role of (i) shadowing, (ii) Cronin effect and (iii) cold nuclear matter energy loss in 
the analysis of open heavy flavor suppression. Our results J^t] bring the treatment of cold nuclear 
matter effects for open heavy flavor production on par with their implementation in the study of 
light hadron and direct photon final states JjJ. 

We now turn to the final-state quark and gluon dynamics. The standard paradigm for jet 
quenching is to consider the energy lost by a parton traversing through the deconfined medium. 
This parton then fragments outside the QGP. We perform the energy loss calculation in the frame- 
work of the GLV approach [7]. The properties of the soft gluon-dominated matter is constrained 
by the experimentally measured charged particle pseudorapidity density ]5|] . 

The results for Raa of 7T° calculated thus, describe the experimental data beautifully. On the 
other hand the experimental results on non-photonic electrons H, ~ R" A ° A for pj > 5 GeV, are 
in clear contradiction with the small quenching of B mesons that give an increasingly important 
(> 50%) contribution to non-photonic e + + e~ in this region (als!]. Cold nuclear matter effects 
increase the discrepancy between the results of partonic level quenching and experiments fl. 

This discrepancy forces a check of our basic assumptions. Whether or not an energetic 
quark traversing the QGP forms a (possibly short-lived) hadronic state in the plasma depends on 
whether the formation time of the hadron, Tf olm jjj, is less or greater than the plasma life time. 
The small pion mass ensures that the parent light quarks and gluons fragment outside of the QGP 
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in accord with the traditional picture of jet quenching. The large D and B meson mass, however, 
implies that charm and beauty quarks will fragment inside the hot medium: Tf 0lm oc 1/m?. Con- 
sequently, the competition between heavy meson dissociation and the c and b quark decay in the 
QGP is a likely physics mechanism that may naturally lead to attenuation of the beauty cross 
section as large as that for charm JUl . 



2. Can a D or B meson exist in the QGP? 

To understand if a D or B meson can exist at rest and in thermal equilibrium in a QGP, and 
to understand how the parton distribution functions (PDFs) and fragmentation functions (FFs) 
are modified in a thermal medium, we solve the Dirac equation for the light quark in a thermally 
modified Cornell potential. Whether such thermal effects do indeed take place for the relevant 
case of an energetic meson moving through the medium, depends on whether the time taken 
to equilibrate this meson is larger or smaller than the time taken to traverse the QGP. We will 
address this question in the next section, and for this section consider the case of a D or B meson 
in equilibrium with the medium. 

Lattice QCD results for the free energy of a system at temperature T containing two infinitely 
heavy quarks separated by a distance r II 1(111 have been used to extract the static potential between 
these quarks. This potential has been employed previously to study the fate of heavy-heavy 
bound states in the QGP formed at temperatures above T c 111 ill . We perform a similar study to 
look at the possibility of the existence of heavy-light mesons in the QGP by looking for bound 
state solutions for the Dirac equation describing the motion of the light quark in the potential set 
up by the heavy quark. The temperature up to which bound states exist depends sensitively upon 
the details of the potential chosen and the effective in-medium mass of the light quark. Here we 
quote results (Table. Q]) for a relatively strongly binding potential and take the in-medium mass 
of the light quark to be the m(T)/ V2. More details are given elsewhere 0. We find that D- and 
B-meson bound-state solutions persist above T c . Their small binding energy and large radius, 
however, will greatly facilitate their subsequent dissolution in the presence of interactions. 

We can use the bound state wavefunctions obtained to calculate the thermally modified PDFs 
and FFs. This calculation is most conveniently done using the wavefunctions in the light cone 



form [12]. The first results for how the PDFs and FFs change with temperature are in ||5|] 



3. Application to heavy meson production in heavy ion collisions 

Heavy flavor dynamics in hot QCD matter critically depends on the time scales involved in 
the underlying reaction. As discussed in the previous section, the dissociation of D and B mesons 
in the vicinity of T c can be facilitated by their small binding energy. Whether such thermal effects 
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Table 1: Properties of the bound D- and B-meson bound-state solutions taking the effective light quark mass to be 
m(T)l V2. These solutions persist up to temperatures ss l.6T c . 
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Figure 1: (Left.) Suppression of D and B hadron production from meson dissociation and heavy quark quenching in 
central Au+Au (top panel) and Cu+Cu collisions (bottom panel) at y.s/vjv = 200 GeV at RHIC. (Right.) Nuclear mod- 
ification for the single non-photonic electrons in central Au+Au and Cu+Cu collisions at RHIC. Data is from PHENIX 
and STAR \2] collaborations. 
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take place in practice, however, depends on the time they need to develop. We can roughly 
estimate this time by boosting the expanded size of the hadron, w 2 yj(r 2 ) from Table [TJ by the y 
factor. When compared to Tf 0Tm this time is large and suggests that the fragmentation component 
of the heavy meson dynamics in heavy ion collisions may not be affected by the QGP. In what 
follows we will study this "instant wavefunction limit". 

A meson that is formed and propagates inside the medium will undergo collisional broaden- 
ing and dissociation [SLL3J]. Following the diagrammatic calculation described in [8[] one can 
calculate the dissociation time Tdi ss for a meson of transverse momentum pj. 

Employing the rates, Tdi ss and Tf orm , the concurrent processes of c and b quark fragmentation 
and D and B meson dissociation are described by rate equations [Ql describing the time evolution 
of quark and meson "densities" (f Q {prJ) = an< i f H (Pr,y,t) = respectively). At 

present, there is no reliable way of incorporating the fluctuations in partonic energy loss in rate or 
transport equations. Therefore, we include the early-time heavy quark inelastic scattering effects 
approximately as a quenched initial condition, f Q (pr, 0) = > f H {Pr, 0) = . Here, the 

attenuated partonic spectrum ^^r~ is calculated differentially versus pj using partonic level 
quenching. The relevant mean quenching time - the time that the physical system of interest 
spends in a quark state - can be estimated from an analytic solution for f e (pr, t), f H (pr, t$\- 
We integrate numerically the above set of coupled ordinary differential equations and use 
the same initial soft gluon rapidity density dN s jdy as in the simulations of 7r° quenching. The 
corresponding suppression of open heavy flavor in y/s^N - 200 GeV central Au+Au and Cu+Cu 
collisions at RHIC is shown in the top and bottom panels on the left in Fig.Q] respectively. For D 
mesons the enhancement from Cronin effect is clearly visible around pj ~ 4 GeV. In both gold 
and copper reactions at RHIC the suppression R B AA as for pj > 4 GeV and these approach 
the quenching of light hadrons for pj > 10 GeV. Future vertex detector upgrades at RHIC will 
ensure direct and separate measurements of the D and B mesons and can be compared directly 
with our results. We also give Raa for non-photonic electrons calculated using semi-leptonic 
decays of charm and beauty hadrons using PYTHIA 11411 . Results for the suppression of open 
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Figure 2: Suppression of D and B meson production in central Pb+Pb collisions at -\Jsnn = 5500 GeV at the LHC 
in two different scenarios. The top panel shows the quenching of heavy hadrons only due to partonic energy loss. The 
bottom panel gives Raa for Ds and Bs with partonic energy loss as well as collisional dissociation of heavy meson. Insert 
shows the corresponding attenuation of non-photonic electrons in a limited pj range. 
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heavy flavor final states in central Pb+Pb collisions at ^snn = 5500 GeV at the LHC taking 
dN s jdy = 2800 are presented in Fig. [2] Raa for non-photonic electrons is shown in the insert. 
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